C-axis Optical Sum Rule in Josephson Coupled Vortex State 
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Observed violations of the c-axis optical sum rule can give important information on deviations 
from in-plane Fermi liquid behavior and on the nature of interlayer coupling between adjacent copper 
oxide planes. Application of a magnetic field perpendicular to these planes is another way to probe 
in-plane dynamics. We find that the optical sum rule is considerably modified in the presence of 
the c-axis magnetic field. Interlayer correlation of pancake vortices is involved in the sum rule 
modification; however, details of the vortex distribution in the plane are less important. 
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I. INTRODUCTION 



The conductivity sum rule is one of the most use- 
ful methods to analyze optical properties of high-T c 
cuprates. The conventional sum ruleEl states that the 
missing spectral weight AN under the real part of the 
optical conductivity between superconducting and nor- 
mal state is equal to the superfluid density p s . However, 
violations of the conventional sum rule along the c-axis 
have been observedB in some high-Tt—cuprates, and are 
related to the change in kinetic energyaQ on entering the 
superconductina-state. Since c-axis response reflects in- 
plane dynamicsJS the c-axis conductivity sum rule and 
the corresponding superfluid density are worthwhile in- 
vestigating to understand characteristics of Cu02 planes 
as well as for their own importance and interest. 

We have previously studiedEl effects of an in-plane mag- 
netic field on the c-axis sum rule and on the super- 
fluid density under the assumption that the in-plane field 
freely penetrates between CuC>2 planes, and found that 
such a field could not easily change the c-axis sum rule. 
The reason for this is that the in-plane field is not di- 
rectly pair breaking and it is only the induced vortex 
due to Josephson tunneling that reduces the superfluid 
density. Usually the conductivity sum rule is not changed 
by a perturbation since the c-axis kinetic energy remains 
unchanged in many cases. In some notion of interlayer 
coupling theory, superconductivity is related to the kk 
netic energy; namely, it might be kinetic energy-driven.B 
In this respect the changes in the kinetic energy and in 
the c-axis sum rule suggest fundamental information on 
superconductivity. Therefore, the c-axis magnetic field 
is more interesting because it gives rise to a shift in the 
quasiparticle energy spectrum, and serves as a direct pair 
breaker. 

When a c-axis magnetic field (H) is above H c \ and 
below H C 2 in a magnitude, where H c i (H C 2) is the lower 
(upper) critical field, and the anisotropy of the Dirac cone 
a d = vf/vq is much larger than one, where vf { v Q^^\ 
the Fermi (gap) velocity, the semiclassical approachEHL2l 
can be applied. Since H ci <C H <C H c2 , vortices are still 
well separated and it is an extended quasiparticle state 
that dominately determines the properties of the in-plane 



dynamics for a d-wave superconductor. For am, ~ 1, 
quantum mechanical effect becomes important ;0 how- 
ever, this is not the case for the high-T c cuprates, for 
which cur) 3> 1. 

In this paper we consider the effect of a magnetic 
field oriented perpendicular to the Cu02 planes on the 
c-axis sum rule and on the corresponding superfluid den- 
sity. We apply a semiclassical approximation to a vor- 
tex state in a two dimensional (i-wave superconductor 
for H c i <C H <C H C 2 and o>d ^ 1, and include in the 
calculation only the energy shift due to the circulating 
supercurrents around the vortex cores. We assume that 
the core occupies only a small part of the single vortex 
unit cell, and that it is the energy shift of the quasiparti- 
cles outside the cores that is most important. We study 
both the the case of correlated and uncorrelated pancake 
vortices, and find that the c-axis sum rule depends on the 
interlayer correlation of pancake vortices while details of 
the intralayer vortex distribution is less significant. 



II. FORMALISM 



We begin with the Hamiltonian H = Ho + 
H c , where Hq describes a d-wave superconductivity, 
and H c .-,= J2i<j *-L [ c ii CT c i2( T + c^Ciia] is interlayer 
couplingjl3 therefore, an electron transfers from a site 
i in plane 1 to the same site i in plane 2. In the semiclas- 
sical approximation, a quasiparticle sitting on the super- 
current circulating around a votex has an energy shift, 
the so-called Doppler shift, and the magnitude of the su- 
percurrent is inversely proportional to the distance from 
the center of the vortex. As illustrated in Figl. (a), for 
the correlated vortex state a quasiparticle sitting on, say, 
the slow supercurrent transfers to the next plane with- 
out change in momentum and, therefore, still sits on the 
slow supercurrent loop as before. Consequently, we ex- 
pect that the kinetic energy will not change during the 
interlayer hopping. Note that the c-axis physics is deter- 
mined by the in-plane dynamics as we pointed out earlier 
[See also Eq. (g)]. On the other hand, for the uncorre- 
lated vortex state, Figl. (b), a quasiparticle on the slow 
supercurrent loop sits on the fast supercurrent after hop- 
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ping; therefore, the kinetic energy changes. The detailed 
calculation requires the in-plane distribution of pancake 
vortices and it will be explained later. 

The optical conductivity follows from the appropriate 
Kubo formula: 



a c (uo) = (i/lu) [lT ret (w) - e 2 d(H c )] 



(1) 



where n ret (w) is the analytic continuation for the 
current-current correlation in Matsubara representation, 
and (H c ) represents the thermal average of the c-axis 
kinetic energy. Here e and d are the electron charge 
and the interlayer spacing, respectively. The super- 
fluid density is detrermined by Kramers-Kronig relation: 
p s = 47rlim w ^o [wlm<7 c (aj)]. In terms of the real part of 
the conductivity o\ c (uj), the c-axis conductivity sum rule 
reads: 



p s = AN — 47re 2 d [(H c ) s - (H c 



(2) 



where AN — 8 J* 0+ c duj [<7i c (a>) — crf c (u;)], and the su- 
perscripts n and s denote the normal and supercon- 
ducting state, respectively. In the integral defining 
AN, the upper limit is sufficiently large so as to in- 
clude all optical transitions of importance to the prob- 
lem but not interband effects. The kinetic energy change 
47re 2 d [(H c ) s - (U c ) n ] will be denoted as AK. Therefore, 
in the presence of a magnetic field (H), the sum rule be- 
comes p s (H) = AN(H) — AK(H), and it can be written 



AN(H) 
Ps(H) 



1 



AK(0) 
P.(0) 



f 



Sp s 

Ps(0) 



SK S 

Ps(0) 



(3) 



where p s (H) = p s (Q) + Sp s and 5K S = AK(H) - AK(0). 
As can be easily seen, the sum rule change due to the 
field depends not only on SK S but also on Sp s . 



III. MAGNETIC FIELD EFFECTS ON THE 
KINETIC ENERGY AND THE SUM RULE 

It is easy to understand that the c-axis magnetic field 
has no effect on the interlayer hopping for a Fermi liq- 
uid in the overdoped regime while the absence of a well- 
established theory of the pseudogap state would seem- 
ingly cause ambiguity in the validity of our approachi 
for the psaudogap regime. However, recently it has been 
observed!!^ that the pseudogap is insensitive to the mag- 
netic field. This implies that even for the underdoped 
cuprates, a semiclassical approach may be applicable to 
the sum rule calculation and the kinetic energy change 
is determined only by the field effect on the supercon- 
ducting state; namely, the Doppler shift. Here, we would 
like to emphasize that we do not calculate the zero field 
kinetic energy difference AK (0) for the pseudogap state 
because a model for the pseudogap would be required 
to do so; insteadj-,we base our analysis on experimen- 
tal observationsntS for this case. We also point out that 



effects of the Doppler shift are considered only on, the 
superconducting state not on the pseudogap stateO 

Let us first calculate the kinetic energy change 5K S . In 
the semicalssical approximation, the fermionic Matsub- 
ara frequencies iu n in the superconducting (2x2) Green 
function G(k, iui n ) are shifted by v s (r) • k, where v s (r) 
is the superfluid velocity at r in a plane and k the quasi- 
particle momentum, that is, the (2 x 2) Green function 
in Nambu space is G(k, ioj n — v s (r) • k), and an average 
over the vortex unit cell is to be carried out. Denoting 
I = v s (r) • k, the (2 x 2) Green function is given by 

A n ■ -a fan - £)t + A k fi + £ k f 3 
G(k, iu n - e) = — — 3^ -2 -2 — • (4) 



fan 



e 2 



Here is a d-wave superconducting gap such that 
the quasiparticle energy = + A 2 , in zero field, 
and fj's the Pauli matrices in spin space. In terms of 
G(k, ioj n — e) , the c-axis kinetic energy in the vortex state 
is given at position r, by K s (r 1 ,r 2 ) = Y^k K s(k; e x , e 2 ), 
where e; = v„i(ri) • k, with 

K s (k;ei,e 2 ) = — ^ Tr |r 3 G(k,MJ- ei)f 3 G(k, iu - e 2 ) 

(5) 

where C — 32ire 2 dt 2 L . The total kinetic energy is obtained 
after K s (r) is appropriately averaged over a vortex unit 
cell of radius R i.e. K S (H) = (1/A) J r<R drK s {r) with 

A = irR 2 the area of the unit cell. Actually, the cor- 
rect averaging procedure depends on whether or not the 
pancake vortices from one plane to the next are corre- 
lated or are completely uncorrelated. For the correlated 
case v s i(ri) = v s2 (r 2 ) and there is a simple space aver- 
age over r while in the uncorrelated case we need to in- 
troduce two separate uncorrelated space variable v s i(ri) 
and v s2 (r 2 ) and take independent averages over ri and 
r 2 i.e. (l/A 2 )/ ri < fl dri/ r3 < fl dr 2 . 

We consider the case of uncorrelated pancakes dis- 
tributed randomly in each plane so that v s i(ri) • k and 
v s2 (r 2 ) • k are to be averaged independently and sepa- 
rately. In a (i-wave superconductor it is the nodal quasi- 
particles that dominate the low T or the small H re- 
sponse. We take T — > and work with the difference in 
kinetic energy with field on and off i.e. K S (H) — K s (0). 
Note that this quantity is determined only by the field 
effect on the superconducting state in our consideration. 
We then apply a nodal approximation. In carrying out 
the sum over k we make a coordinate transformation to 
p with the nodal point k„ on the Fermi surface play- 
ing the role of the origin for p. We then integrate over 
a nodal region and the integration over p should not be 
sensitive to the region of integration provided |p| is taken 
to be of order po ~ O(Ao). With these simplification we 
can introduce a (5-function S (ei — v s i(ri) • k„) into the 
expression for K S {H) — K s (0) and define 
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V{ex) = — J dri5 (e x - v s i(ri) • k„) , 



(6) 



to obtain 

/oc 
de 1 de 2 V(e 1 )V(e 2 ) 
-OC 

x[*(ei,e 2 )-*(0,0)] , (7) 

where ^(ei^) = 2 k if 8 (k;ei,e2). The matrix sum in 
^(ei, 62) can be done to get 



C, 



G(k, iuj n - ei)G(k, iw n - e 2 ) 



-F(k, iw„ - ei)F(k, iuj n - e 2 ) 



(8) 



with G(k, zw„) and -F(k, iw„) the superconducting state 
ordinary and anomalous amplitude, respectively. For 
simplicity we will only treat the pure case, and after some 
algebra we derive ^(ei, e 2 ) as follows: 




(9) 



where fi± — /(ej ± Eu) is the Fermi function, and 
(12 = ei — e 2 - Note the symmetry properties: ^(ei, e 2 ) = 
*(e 2 ,ei) = *(-ei,-e 2 ) and *(-ei,e 2 ) = *(ei,-e 2 ). 

The limit ei — > e 2 is of interest because it enters the 
correlated pancake case for which v s i(ri) k = v s2 (r 2 )-k. 
In this case the average for the kinetic energy over a sin- 
gle V(e) is 

/oo 
d£P(e)[#(e,e)-¥(0,0)] , (10) 
-00 



where 



C 



f2 



4 4^ F 2 



02 



df(z) 



dz 



/(e + £ k )-/(e-£ k ) 



(11) 



In the nodal coordinate system E'k = \J p\ + P 2 = P an d 
^ k = X^node ^ / pdpdfi, where 1? is a polar angle and 

J = [(2tt) 2 vfvg\ ■ At zero temperature the Fermi 
function become 0-function i.e. f[z) — > 0(— 2) such 
that 0(— 2) — for 2 > and otherwise G(— 2) = 1. 
Substututing ^(e, e) into Eq. jl0|), we obtain 

K s (H)-K s (0) = 2C deV(e)J J pdpdd 



Pi 



-5(e+p)-5(e-p) 



9(-e-p)-0(-e + p) + l 



(12) 



The first (5-function and 0-function in Eq. ( |l2| ) do not 
survive as they give zero contribution. The first integral 
in Eq. (|T^) over p gives a contribution proportional to e. 
The second integral is limited to the range pe (e,Po), but 
its range can be extended to (0,po) by including the con- 
tribution of the first integral. This combined contribu- 
tion cancels the last integral in Eq. (|l^) to give zero. As 
pictorially illustrated in Figl. (a), we analytically showed 
that K S {H) = K s (0) i.e. the out- of -plane magnetic field 
has no effect on the c-axis kinetic energy regardless of 
the form ofV(e) in the correlated vortex state. This does 
not mean, however, that there is no change in the sum 
rule given by Eq. (^). Instead it means that the sum 
rule is modified only by a possible change in superfluid 
density. The calculation of 5p s is similar to the above 
[See later explanation] . It turns out that p s (H) — p s (0) = 
AttCJ J °° de eP(e), where p s (0) = CN(0)/2 with the den- 
sity of states at the Fermi level N(0), and we obtain 

Wp*(oh-(2/a ) ;;°cj, e p( e ). 

Recently, Vekhter et.al£3 have given analytic expres- 
sions for different models of the vortex distribution func- 
tion Vie). They consider a Gaussian Vg{c) and two pos- 
sible vortex liquid model 'pLi(e) and 7- > i 2 (e): 



V L x{e) 
V L 2(e) 



1 



1 



exp 



El 



2{E 2 H + e 2 ) W 2 ) 
1 



(13) 
(14) 



ttE h 



y/(2e/E H ) 2 + 1 



e 3 4e 5 



3E% 
2e 4 



(15) 



The magnetic energy Eh is given by vf/(2R), where 
the radius of the vortex unit cell R = -J 1 Qq/ttH with 
$0 the flux quantum. Using the model V(e), we obtain 
6p s /p s (0) = —aiEjj/Ao. where a.\ is for G, 1 for 

LI, and 0.84 for L2. Note that all distributions give the 
same linear dependence on Eh / Ao with slightly different 
coefficients, and H reduces the superfluid density as we 
expect. 

Our main results so far are that in the correlated vortex 
case there is no change in kinetic energy on application 
of a magnetic field H perpendicular to the layers, inde- 
pendent of the vortex distribution function V(e). In all 
cases considered, the superfluid density is reduced by the 
magnetic field, and Sp s /p s (0) is linearly proportional to 
Eh/Ao with a coefficient of order 1. The form of V(e) 
determines the exact value of the coefficient. We con- 
sider two interesting regimes of doping for the correlated 
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vortex case (SK S = 0). For the overdoped case we as- 
sumed that the in-plane motion can be described by a. 
Fermi liquid. In this case, it has been found previouslyO 
that in zero field AiV(0)/p s (0) = 1 i.e. the c-axis 
sum rule is conventional, which is consistent with ex- 
perimental observations!] It arises because AK(0) = 
in this instance. Therefore, it follows from Eq. (||) that 
AN(H)/ p s (H) — 1 also i.e. there is no change in the 
sum rule induced by the external field. 

As we mentioned earlier, we do not rely on a spe- 
cific theory to calculate the sum rule for the under^. 
dope regime. Instead we use experimental observationscl 
on the underdoped cuprates. In this case the zero 
field sum rule is observed to be about a half; namely, 
AA(0)/p s (0) ~ 1/2, which also follows directly in the 
preformed pair model for the pseudogap. This means 
that AK(0)/p s (0) = -1/2 in Eq. @. When we calculate 
SK S , we take into account only effects of the Doppler shift 
due to the—magnetic field, again, based on experimen- 
tal resultsES Consequently, we obtain AN{H)/ p s {H) = 
1/2 + Sp s /2p s {0) = 1/2 - ai E H /2A , where a x = 
for G, 1 for LI, and 0.84 for L2. The sum rule originally 
equal to 1/2 is further reduced by the magnetic field. 

Now, we return to Eq. (|?]) for the uncorrelated vor- 
tex state with Eq. (^|) defining 3>(ei,e 2 ). In this case 
K s (H)~K s (0) becomes 

K s (H) - K s {0) 

= 2 / de 1 de 2 V(e 1 )V(e 2 ) [*( £l , e 2 ) - (0, 0)] 
Jo 

C °° 

+2/ deide 3 P(ei)P(e 3 )[¥(e 1 ,-e a )-¥(0,0)] . (16) 
Jo 

On application of the nodal approximation to the differ- 
ence \?(ei, e 2 ) — W(0, 0) in the limit T — *■ 0, we obtain 



6K a = 



CN(0) 



2 A 

x 



/>0O 

/ de 1 de 2 V{e 1 )P(e 2 ) 
Jo 



6162 , l ( i \ 1 

ei+£2) 2 



ei + £2 



(17) 



which is of order Eh/Aq, and numerical calcu- 
lations give SK s /p s (0) = — a 2 En /A , where a 2 
is 0.123 for G, 0.135 for LI, and 0.102 for L2. 
We also derive the change in superfluid density 
which follows in a similar fashion; namely, 5p s — 
de x de 2 V{e l )V{e 2 ) [ X {ei, e 2 ) - x(0, 0)], where 



X(£i,£2) 



C 



Al 



fl+ - f-2+ , fl- - /2 



4 ^ El 



fi+- h 



£12 



fl- - h 



£12 



ei2 + 2^ £12 — 2i?k 



(18) 



It can be shown that the superfluid density reduction 
becomes 



Sp s 



C N(0) 
'2 A 



de 1 de 2 P(e 1 )7'(e 2 ) 



3(ei + £ 2 ) 2 + ej + e| (ei + e 2 ) 

2(£l+£ 2 ) 2 



In 



£i - £2 



£i + £2 



and Sp s /p s (0) = —asEn/Ao, where is 0.34 for G, 
0.72 for LI, and 0.63 for L2. For correlated vortex, we 
used X (£,£) with x(0,0) = p s (0). 

The sum rule of the uncorrelated vortex state reads 
AN(H)/ p s (H) ~ l-a 2 (E H /A ) for the overdoped, and 
AN(H)/p s (H) ~ 1/2- (a 2 + a 3 /2)(£; H /A ) for the un- 
derdoped regime. In both cases the sum rule is reduced 
by the presence of a magnetic field. In the overdoped case 
the reference is, however, 1 while for the underdoped case 
it is 1/2, and the coefficient of the reduction, which goes 
like Eh/Aq, is larger for the underdoped case as com- 
pared with the overdoped case. Note that the magnetic 
energy Eh goes like \fll in all cases. 



IV. CONCLUSIONS 

For correlated pancake vortices with an in-plane Fermi 
liquid and coherent c-axis coupling, there is no change 
in the c-axis optical sum rule when a magnetic field is 
applied perpendicular to the layers of a d-wave supercon- 
ductor. The sum rule keeps its conventional value of one. 
At the same time the superfluid density in the c-direction 
decreases by an amount proportional to the square root 
of the magnitude of the external magnetic field (y/~H). If, 
however, the in-plane dynamics is unconventional, for ex- 
ample, it can be described by a pseudogap above T c , the 
sum rule, which is equal to 1/2 when H = 0, gets reduced 
linearly in y/~H as does the c-axis superfluid density. For 
uncorrelated vortices, there is always a reduction in the 
sum rule which is also proportional to \fli . For the in- 
plane Fermi liquid case the sum rule is reduced below 
one, while for the unconventional case with a pseudogap 
it is reduced below 1/2. In both cases a reduction in 
the c-axis superfluid density accompanies the sum rule 
reduction. 
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FIG. 1. Pictorial illustration for a correlated (a) and an uncorrelated (b) vortex state. v s is the supervelocity and the gray 
tubes are vortex cores. 
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(a) Correlated vortex (b) Uncorrected vortex 



